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SUMMARY 


The  objective  was  to  develop  an  algorithm  for  describing 
elastic-viscoplastic  constitutive  equations  based  on  the  Bodner- 
Partom  model  in  the  STEALTH  finite  difference  code.  To  ensure 
convergence  of  the  B-P  solutions,  several  parameters  (TOLRSD, 
RLAXSD,  MITRSD,  MCUTSD,  and  TCFSD)  were  introduced  into  the 
algorithm.  Based  on  the  results  of  various  test  cases,  the 
algorithm  seems  to  work  efficiently  using  the  default  values  of 
these  parameters.  However,  when  convergence  problems  arise,  the 
user  can  adjust  the  values  accordingly. 

The  validity  of  the  newly  developed  routines  (described  in 
Appendix  A)  were  checked  through  two  different  methods.  In  the 
first  method,  a  dynamic  tensile  test  under  uniaxial  stress  state 
was  simulated  using  STEALTH.  The  stress-strain  responses  at 
various  strain  rates  were  successfully  compared  with  the  direct 
solutions  of  the  Bodner-Partom  equations  under  uniaxial  stress 
state.  In  the  second  approach,  for  a  rate- independent  elastic- 
perfectly  plastic  material,  a  plate  impact  test  simulation  was 
carried  out  using  (1)  the  regular  material  model  that  is 
available  in  STEALTH,  and  (2)  the  developed  Bodner-Partom  model. 
The  results  in  terms  of  velocity  and  stress  history  in  the  target 
compared  extremely  well. 

We  have  used  quite  successfully  the  Bodner-Partom  model 
option  that  has  been  described  in  this  report  to  characterize 
various  metals.  We  can  describe  rate  dependent,  rate- 
independent,  and  strain  hardening  material  behaviors  using 
STEALTH  with  the  new  material  option;  the  corresponding  results 
for  different  steel,  aluminum  and  copper  are  given  in  Reference 
10.  The  development  and  the  demonstrated  use  of  a  sophisticated 
strain-rate  and  history  dependent  material  model  in  a  finite 
difference  code  will  certainly  encourage  many  code  users  who  work 
in  armor  design  to  employ  such  a  model  in  the  computer  simulation 
of  high  velocity  impact  problems. 
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SECTION  1 
INTRODUCTION 

Several  applications  of  structural  materials  involve  impact 
type  loading.  Under  such  loading  conditions,  the  material 
experiences  high  rates  of  strain,  large  deformation,  high 
pressures  and  high  temperatures.  The  dynamic  materials 
characterization  becomes  very  important  in  the  analysis  of  weapon 
effects,  including  structural  response  to  explosive  loading. 

Within  the  last  decade,  the  development  of  finite 
element/difference  codes  has  provided  additional  analytical 
capability  in  understanding  these  problems.  It  is  now  well 
established  (see  Reference  1)  that  the  material  descriptions  can 
significantly  affect  the  computational  results.  Code  users  have 
often  been  unable  to  apply  sophisticated  constitutive  relations 
for  describing  materials.  For  example,  most  computer  codes  still 
do  not  account  for  the  strain-rate  and  pressure  dependence  of 
yield  and  flow  stresses. 

Several  constitutive  theories  based  on  state  variables  have 
been  developed  by  various  investigators  (References  2  thru  6).  One 
of  the  advantages  of  state  variable  theories  is  that  changes  in 
constitutive  behavior  with  load  history  can  be  modelled  by 
suitably  modifying  the  evolution  equation  of  the  corresponding 
state  variable.  Depending  on  the  complexities  of  the  loading, 
more  than  one  state  variable  may  be  introduced  into  the  theory. 

The  objective  of  the  present  work  is  to  develop  algorithms 
(computer  subroutines)  describing  rate-dependent  (viscoplastic) 
constitutive  equations  for  a  general  purpose  finite  difference 
code.  For  this  purpose,  a  state-of-the-art  finite  difference  code 
called  'STEALTH'  (see  Reference  7)  was  selected. 

The  state  variable  based  viscoplastic  theory  developed  by 
Bodner  and  Partom  (see  Reference  2)  was  chosen  in  the  present  work 
for  three  reasons:  (1)  ability  to  predict  the  response  of  material 
to  a  broad  range  of  load  histories,  (2)  adaptability  to  finite 
difference  analyses  of  structural  components,  and  (3)  ease  with 
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which  the  parameters  in  a  constitutive  model  can  be  determined 
from  high  strain-rate  experiments. 

the  Bodner-Partom  model  is  described  in  Section  2.  The 
yield  models  in  STEALTH  and  the  new  B-P  model  algorithms  are 
briefly  discussed  in  Section  3.  The  validation  of  the  developed 
routines  is  discussed  in  Section  4.  Finally,  Appendix  A  discusses 
the  additions  required  to  include  the  Bodner-Partom  model  in 
STEALTH,  while  Appendices  B,  C,  and  D  supply  additional  details 
concerning  the  input  preparations,  new  subroutines,  and  flow¬ 
charts.  To  understand  the  modifications,  it  is  essential  that  one 
is  familiar  with  the  STEALTH  code  structure.  We  have  strictly 
followed  the  STEALTH-code  conventions  in  the  input  preparations. 
Since  the  STEALTH  manual  (Reference  7)  describes  in  detail  the 
conventions  that  have  been  employed  in  its  preparations,  we  have 
described  only  the  additions  that  are  needed  to  use  the  Bodner- 
Partom  model.  The  iterative  solution  scheme  that  has  been  used  in 
the  Bodner-Partom  model  routines,  in  general,  can  be  understood  by 
any  code  develop er/user  from  the  self-explanatory  flowcharts  in 
Appendix  D. 
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SECTION  2 

BODNER-PARTOM  MODEL 

For  assumed  small  strains,  the  total  strain  rate  is  taken  to 
be  decomposable  into  elastic  (reversible)  and  inelastic 
(nonreversible)  components 


e .  . 

ij 


=  e 


1J 


*P 

+  e*:  . 
ij 


(1) 


which  are  both  nonzero  for  all  loading/unloading  conditions.  The 


elastic  strain  rate,  e.  ., 

^  J 


is  related  to  the  stress  rate  by  the 
time  derivative  of  Hooke's  Law  (based  on  the  assumption  of  small 
strains).  The  inelastic  strain  rate  can  be  expressed  in  the 

general  form 


•p 

e*: 


ij 


ePjtCij.  Zk>T) 


(2) 


where  Zk  are  one  or  more  internal  (inelastic)  state  variables  and 
T  is  the  temperature.  In  particular,  e1?  .  is  taken  to  follow  the 
Prandtl-Reuss  flow  law  of  classical  plasticity,  (assuming 
incompressibility) : 


e?  .  =  e?  .  =  AS.  . 
ij  ij  ij 


(3) 


•  p 

where  et^  and  are  the  deviatoric  plastic  strain  rate  and  the 

stress,  respectively.  Squaring  equation  3  gives 

1/2  e^  e*?j  =  =  1/2A2S.j  S.  ^  =  A2J£ 

or 


dP  =  a2j2 


(*0 


where  is  the  second  invariant  of  the  plastic  strain  rate  and  J2 
is  the  second  invariant  of  the  stress  deviator. 
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The  basic  assumption  of  the  formulation  is  that  all 
inelastic  deformations  are  governed  by  a  continuous  relation 
between  and  J2.  The  particular  form  that  was  established  to 
describe  this  relation  was  based  on  the  concepts  and  equations  of 
"dislocation  dynamics"  and  it  is  given  by: 


{exp  [(- 


Z2  Nn  (n  +  1)  i ! 
3J2;  n  JI 


(5) 


where  Dq  is  the  limiting  value  of  the  plastic  strain  rate  in 
shear.  n  is  a  material  constant  that  controls  strain  rate 
sensitivity  and  also  influences  the  overall  level  of  the  stress- 
strain  curves.  It  is  a  fundamental  constant  and  it  is  not 
dependent  on  the  loading  history.  Z  is  the  inelastic  state 
variable  which  represents  the  measure  of  the  overall  resistance  to 
plastic  flow  caused  by  microstructural  barriers  that  impede 
dislocation  motion.  Combining  equations  3  through  5,  the  plastic 
strain  rate  takes  the  form 


Dq  exp 


n 


(n+1)  i 
2n  1 


(6) 


The  general  form  for  the  evolution  equations,  i.e., 
dependence,  of  the  inelastic  state  variable  Z  is, 


history 


Z  =  F(J2,  Z,  T)  (7) 

For  conditions  of  uniaxial  stress  of  constant  sign,  the 

hardened  state  with  respect  to  plastic  flow  is  assumed  to  be 

represented  by  a  single  state  variable  Z  which  depends  on  plastic 

work.  This  corresponds  to  isotropic  hardening  and  the  evolution 

equations  employed  can  be  shown  to  be  in  the  form  of  equation  7. 

The  corresponding  form  is  then  based  on  the  concept  that  plastic 

work,  Wp ,  controls  the  hardening  process,  and  that  the  plastic 

work  and  its  time  derivative  are  functions  of  o . .  and  e. ..  The 

i  J  i  J 
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evolution  equation  for  the  state  variable  is  assumed  to  have  the 
form : 

Z  =  m(Z1-Z)  Wp  (8) 

where,  W  «*  o.  .  ep..  The  above  equation  can  be  integrated  With  m 

r  A  J  A 

and  Z,  as  material  constants.  The  relation  between  Z  and  Wp 
becomes, 

-mW 

Z  =  Z1-(Z1-ZQ)  e  P  (9) 

ZQ  is  the  initial  value  of  Z  corresponding  to  the  reference  state 
from  which  Wp  is  measured.  It  is  noted  that  ZQ  could  take  any 
value  between  0  and  Z^.  The  initial  yield  stress  of  the  uniaxial 
stress-strain  relation  would  be  directly  influenced  by  ZQ.  Z-^ 
corresponds  to  the  saturation  state  of  Z  at  which  the  material 
reaches  a  nonwork  hardened  condition,  m  is  a  material  constant 
that  controls  the  rate  of  work  hardening. 

For  materials  which  exhibit  strong  strain  hardening,  m  was 
made  a  function  of  Wp  by  Bodner  and  Merzer  (Reference  9).  The 
expression  for  m  was: 

-aW 

m=m+m^ep  (10) 

oi 

which  adds  two  more  constants.  When  m  is  defined  by  equation  10, 
the  relationship  (equation  9)  is  not  valid.  The  expressions  for  Z 
when  integrated  after  replacing  m  by  the  definition  (equation  10), 
yields  the  following  relationship  between  Z  and  Wp : 

-m  W  -(m  +m1-m) 

Z  -  Z1-(Z1-Z0)  e  P  e  -  (11) 
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•  p 

The  plastic  strain  rate  expression,  (equation  6)  e^j, 
corresponding  to  uniaxial  stress  state  ( ox  4  0,  °y  =  °z  =  ^)» 
becomes 


D-  exp  {-  ( 


n+1 

2n 


7  2n 

>(|)  1 


(12) 


The  expression  for  the  state  variable,  Z,  becomes 


-mW 

z  =  Z1  -  (Z^Zq)  e  P  (13) 

For  a  given  stress  or  strain  rate  history  under  uniaxial  stress 
state,  the  calculation  of  stress/strain  response  is  fairly 
straight  forward  through  a  first  order  Euler's  integration  scheme 
of  equations  (12)  and  (13). 

In  general,  the  stress/strain  response  under  a  multiaxial 
stress  state  requires  the  solution  of  the  second  order  tensor 
as  expressed  by  the  equation  (6).  The  algorithm  developed  for  the 
STEALTH  code  is  for  three  dimensional  states  of  stress.  The 
essential  changes  that  were  introduced  in  the  regular  STEALTH 
routines  are  discussed  in  detail  in  Appendices  A  through  D.  The 
following  section  discusses  the  results  obtained  from  STEALTH  for 
the  validation  of  B-P  model  algorithms. 
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SECTION  3 

INCORPORATION  OF  B-P  MODEL  ALGORITHM 

The  elastic-plastic  constitutive  model  in  the  standard 
STEALTH  version  is  not  capable  of  describing  time  dependent 
material  response.  The  description  of  the  yield  model  in  STEALTH 
is  fairly  straight  forward.  In  this  section,  we  briefly  discuss 
the  algorithm  that  is  available  in  STEALTH  to  describe  rate 
independent  elastic-plastic  material  response.  The  B-P  model 
algorithm  to  describe  strain  rate  and  history  dependent  material 
response  is  also  discussed. 

3.1  YIELD  MODELS  IN  STEALTH 

Several  yield  models  in  which  the  flow  stress  is  function  of 
pressure,  internal  energy,  and  distorsional  energy  are  available 
in  STEALTH.  Those  models  are  in  general  empirical  and  they  can 
not  describe  the  strain  rate  and  loading  history  dependent  yield 
and  flow  stresses.  The  simplest  model  that  is  available  to 
calculate  the  flow  stresses  which  do  not  depend  on  pressure  is  the 
elastic-perf ectly  plastic  material  model.  The  model  is  given  by, 

3J2  =  Y2 

where  Y  is  the  constant  yield  stress.  Y  is  the  only  material 
parameter  that  is  needed  to  describe  the  model. 

In  STEALTH,  the  stress  deviators,  are  estimated  based  on 

the  elastic  material  behavior  at  the  end  of  each  time  step.  The 
effective  stress,  is  estimated  from  the  following 

relationship : 


o 


eff 


°eff  is  then  compared  with  the  yield  stress  (Y).  When  oeff  is 
greater  than  Y,  the  plastic  flow  is  assumed  to  occur.  The  values 

0  g  £ 

of  the  estimated  are  multiplied  by  a  correction  factor  (3)  so 
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that  the  S.  .  lie  on  the  Von-Mises  yield  surface.  The 

J 

corresponding  expression  is  given  by. 


3  can  be  calculated  as: 


The  stress  deviator  S^  ^  corresponding  to  the  plastic  flow  is  then 
given  by, 


sij  - 


When  there  is  no  plastic  flow,  i.e.,  3  is  equal  to  1.  The 

important  assumption  in  the  above  algorithm  is  that  the  stress 
deviators  do  not  depend  on  the  loading  rate,  time  or  history.  As 
mentioned  earlier,  when  the  material  behavior  is  viscoplastic, 
these  kinds  of  yield  models  are  no  longer  valid.  In  fact,  most  of 
the  metals  exhibit  rate  dependency  under  very  high  strain  rates. 
The  calculations  of  S .  .  will  require  an  iterative  scheme.  Since 

A  J 

the  B-P  model  is  a  viscoplastic  model,  we  developed  an  algorithm 
in  which  the  stress  deviators  are  computed  through  a  subincre- 
mental  time  stepping  scheme  as  described  in  the  following  section. 


3.2  BODNER-PARTOM  MODEL  IN  STEALTH 

We  described  the  B-P  model  earlier  in  detail  (Section  2). 
The  B-P  model  is  simply  given  by  the  following  two  relationships: 


D 


C  n  +  1 
l2n 


n 

) 


(14) 
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and 


Z  =  m(Z1-Z)W 


(15) 


A  sub  incremental  time  stepping  scheme  is  used  to  calculate 
the  deviatoric  stresses,  for  each  finite  difference  time-step. 

The  global  time  step.  At,  is  further  divided  into  small  steps  in  a 
special  purpose  subroutine  to  describe  the  B-P  model.  At  the  end 
of  each  global  time  step,  STEALTH  provides  the  total  strain  rates, 
e^,  the  volumetric  strain  rate,  e,  and  the  pressure,  P;  S^,  °ij» 
and  are  unknown. 

In  the  current  numerical  scheme,  the  plastic  strain  rates 
are  estimated  for  each  sub  incremental  time  step  using  equations  14 

and  15  with  the  value  for  J?  calculated  from  the  known  S .  ,  values 

•  6 

at  the  beginning  of  the  step.  The  elastic  strain  rates,  e..,  are 

*  D  ^ 

estimated  from  the  estimated  and  the  known  total  strain  rates, 

e.  ..  The  new  estimates  for  S,  .  can  be  made  from  the  following 
ij  iJ 

relationship : 


S .  ,  =  o..  +  p5  .  . 
ij  ij  y  1J 


(16) 


where 


o  .  .  *  E . ,  e!f  . 
ij  il  lj 


(17) 


are  the  elements  of  the  elastic  modulus  matrix. 

Using  the  new  estimates  for  S^,  the  improved  values  for  the 
plastic  strain  rates  at  the  end  of  the  sub  incremental  time  step 
can  be  computed.  This  procedure  is  continued  until  the  values  of 
5 ^ j  computed  for  the  two  successive  iterations  converge  to  the 
same  values  within  the  tolerant  limits. 
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The  algorithm  has  been  designed  to  handle  convergence 
problems.  If  any  one  component  of  the  stress  deviators,  S.  .,  does 
not  converge  within  the  imposed  iteration  limit  and  convergence 
tolerance,  the  sub  incremental  time  step  is  reduced.  The  solution 
procedure  for  the  current  time  interval  is  repeated  using  this 
smaller  time  step.  If  the  still  fail  to  converge  after  the 

maximum  allowable  number  of  time  step  reductions,  the  program 
creates  a  debug  output  file  for  further  diagnosis  and  stops 
execution. 

The  user  can  also  control  several  other  parameters  in  the 

convergence  of  the  Bodner-Partom  solution.  A  relaxation  factor,  R 
% 

(RLAXSD)  ,  is  used  to  weigh  the  previous  estimates  with  the 
current  in  order  to  obtain  improved  values  for  the  next 

iteration,  as  follows: 

Sest  =  R • S .  .  +  (1-R)*S®®fc  (18) 

1J  ij  1J 

When  R=1 ,  the  new  estimates  are  equal  to  the  current  values  and 

when  R=0 . 5 ,  the  new  estimates  are  simply  the  averages  of  the 

previous  estimates  and  the  new  values.  Generally,  R=1  works  well; 

sometimes,  however,  R  must  be  reduced  (to  0.75  or  0.5)  in  order  to 

prevent  an  oscillating  solution.  Another  parameter  which  affects 

£ 

convergence  is  Ae  (DEMAX)  ,  the  maximum  allowable  change  in  the 

total  strain  during  a  sub  incremental  time  step.  A  reduction  in 

Aemax  can  result  in  a  smaller,  initial  sub  incremental  time  step, 

which  should  help  the  S.  .  converge  more  rapidly.  We  can  also 

1 J  * 

adjust  the  relative  convergence  tolerance  (TOLRSD)  to  help 

achieve  convergence.  Likewise,  an  increase  in  the  maximum  number 

* 

of  iterations  (MITRSD)  allowed  may  sometimes  give  the  solution  a 

better  chance  of  converging.  Finally,  a  decrease  in  the 

sub inc rementa 1  time  step  size  may  aid  in  the  convergence  of  the 

S.  . ;  this  can  be  accomplished  by  increasing  either  the  maximum 
1 J  * 

number  of  time  step  reductions  (MCUTSD)  or  the  time  step 

* 

reduction  factor  (TCFSD)  .  These  various  convergence  related 


*  . ~ " . . 

See  Appendix  C 
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parameters  are  described  in  Appendix  C  as  input  parameters.  We 
found  that  several  test  problems  converged  consistently  using  the 
default  values  of  these  parameters.  A  user  may  use  the  default 
values  unless  convergence  problems  are  encountered. 

To  incorporate  the  B-P  model  into  STEALTH,  we  developed 
several  subroutines.  A  series  of  flowcharts,  given  in  Appendix  D, 
describes  the  steps  involved  in  the  numerical  scheme.  The 
following  section  focuses  on  the  model  problems  which  were 
considered  for  the  validation  of  the  developed  numerical  scheme. 
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SECTION  4 

VALIDATION  OF  THE  ALGORITHM 

The  new  subroutines  that  describe  the  Bodner-Partom 
constitutive  equations  must  be  checked  for  their  validity.  For 
this  purpose,  two  approaches  were  considered.  In  the  first 
(direct)  approach,  the  stress-strain  response  of  a  uniaxial 
tensile  bar  was  computed  from  STEALTH  using  the  new  routines.  The 
results  were  compared  with  the  solution  obtained  from  an 
independent  program  (BPSOLVE)  which  solves  the  uniaxial  Bodner- 
Partom  equations  (12)  and  (13)  using  Euler  forward  integration 
scheme.  Thd  second  approach  is  an  indirect  way  of  checking  the 
validity  through  a  plate-impact  test  simulation  (see  Reference  8). 
In  this  section,  these  two  approaches  are  discussed  in  detail. 

4.1  DIRECT  APPROACH 

A  thin  cylindrical  bar  being  pulled  in  tension  at  a  constant 
velocity,  as  shown  in  Figure  1,  was  considered.  The  bar  was 
modeled  using  the  'THIN  OPTION'  that  is  available  in  STEALTH  (see 
Reference  7).  This  option  provides  solutions  to  problems  in  which 
the  stress  state  is  one-dimensional  as  in  the  uniaxial  tensile 
test.  In  the  finite  difference  analysis,  the  bar  was  modeled  with 
10  zones. 

We  considered  two  different  material  behaviors:  (1)  rate 
independent,  elast ic-perf ectly  plastic  solid,  and  (2)  rate 
dependent,  elastic-viscoplastic,  strain  hardening  solid.  The  B-P 
model  parameters  used  in  the  first  case  were: 

Dq  =  108/sec 
ZQ  =  15  kbar 
Z1  =  15  kbar 
m  =  0.2  kbar  1 
n  -  6 . 5 

Typically  a  large  value  of  n  indicates  that  the  material  is  rate 
independent.  The  same  initial  and  saturation  values  of  Z^ 
indicate  that  the  material  is  perfectly  plastic.  In  this 
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particular  case,  the  strain  hardening  parameter  m  cannot  influence 
the  material  response  and  it  becomes  a  redundant  parameter. 
However,  an  arbitrary  value  of  0.2  was  used  for  m  in  the 
simulation.  Using  the  Bodner-Partom  model  algorithm  in  STEALTH, 
the  stress  response  in  the  tensile  bar  was  computed.  The 
corresponding  input  data  for  the  STEALTH  run  is  given  in  Appendix 
E.  The  pull  velocity  in  the  run  was  assumed  to  be  10  m/sec.  The 
average  strain  rate  in  the  bar  corresponding  to  this  velocity  is 
100/sec.  This  average  strain  rate  is  typically  experienced  by 
zone  5  which  is  the  middle  zone.  The  stress-strain  curve  obtained 
from  BPSOLVE  for  a  strain  rate  of  100/s  compared  extremely  well 
with  the  STEALTH  results  as  can  be  seen  from  Figure  2.  We  also 
calculated  stress-strain  curves  for  strain  rates  1000  and 

_  i 

10,000s  .  The  curves  were  essentially  the  same  as  in  Figure  2. 

This  validates  the  solution  algorithm  since  the  material  was 
modeled  as  a  rate  independent  solid. 

We  repeated  the  simulation  of  the  uniaxial  tensile  bar  test 
modeling  the  material  as  a  rate  dependent,  strain  hardening, 
elastic-viscoplastic  solid.  The  model  parameters  we  chose  for 
this  case  were: 

D  =  10^/sec 
o 

Z  =  55  kbar 
o 

Zn  =  70  kbar 

1  -1 
m  =  1.5  kbar 

n  =  0.4 

A  relatively  low  value  of  n  indicates  the  response  is 

significantly  rate  dependent.  The  strain  hardening  behavior  is 

represented  by  the  different  values  for  ZQ  and  ,  and  a  value  of 

1.5  for  m.  The  cylindrical  bar  problem  was  simulated  using 

STEALTH  with  these  B-P  model  parameters.  Three  different  runs 

were  made  with  pull  velocities  of  10,  100,  and  1000  m/sec.  The 

2  3 

corresponding  average  strain  rates  in  the  bar  were  10  ,  10  ,  and 
4 

10  /sec.  Using  BPSOLVE  the  stress-strain  responses  for  these 
strain  rates  were  obtained  and  were  compared  with  the  STEALTH 


•H  .«  .t*  .IS  .30  .35 

•TRAIN 


Figure  2.  Comparison  of  Stress-Strain  Curves  Obtained 
from  STEALTH  Using  B-P  Routines  and  BPSOLVE 
at  a  Strain  Rate  of  100  s“l  for  Rate-Independent 
Material  Behavior. 
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results  as  shown  in  Figure  3-  It  can  be  seen  that  the  STEALTH 
solutions  [solid  line]  compared  extremely  well  with  the  direct 
solutions  [symbols]  of  the  B-P  equations. 

4.2  INDIRECT  APPROACH 

The  idea  behind  this  approach  was  to  compare, the 
stress/velocity  history  of  the  target  plate  obtained  from  a 
conventional  plate-impact  test  simulation  (see  Reference  8  for 
details)  using  two  constitutive  models:  (1)  the  rate- independent 
elast ic-perfectly  plastic  material  model  that  is  already  available 
in  STEALTH,  and  (2)  the  Bodner-Partom  model  that  is  included  in 
STEALTH  as  a  special  purpose  material  model.  The  rate  independent 
elastic-perf ectly  plastic  model  used  in  the  first  simulation  was 
modeled  by  Bodner-Partom  equations  and  the  corresponding  B-P  model 
constants  were  used  in  a  second  simulation  of  the  same  problem. 
Since  the  material  behavior  used  in  both  simulations  is 
essentially  the  same,  the  velocity  and  stress  history  of  the 
target  calculated  from  the  STEALTH  with  and  without  the  B-P  model 
should  be  the  same.  This  indirect  approach  was  used,  and  results 
were  obtained  for  an  impact  velocity  of  300  m/sec.  The  input  data 
used  in  STEALTH  is  given  in  Appendix  F. 

The  B-P  model  constants , corresponding  to  the  rate- 
independent  elastic-perf ectly  plastic  materials  were, 

Zq  =  15.0  kbar 
Z.  =15.0  kbar 

1  -i 

m  =  0.2  kbar 
n  =  6.5 
Dq  =  1 08/sec 

The  uniaxial  stress-strain  curve  that  is  described  by  the 
B-P  model  using  these  constants  is  shown  in  Figure  4.  The  flow 
stress,  Yq,  in  the  figure  was  employed  in  the  regular  STEALTH  run 
(without  B-P  model  routines)  and  the  corresponding  value  was  Yq  = 
11.8  kbar. 
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STRESS 


The  Elastic-Perf ectly  Plastic  Stress-Strain 
Curve  [ —  Solid  Line]  and  B-P  Simulation 
[Symbol]  used  in  STEALTH  Simulations  of 
Plate  Impact  Tests. 


Figure  4 . 


As  can  be  seen  from  Figures  5  and  6,  the  free  surface 
velocity  of  the  target  and  the  stress  history  at  the  spall  plane 
obtained  from  the  two  cases  are  almost  identical.  This  further 
validates  the  newly  developed  subroutines  that  describe  the 
Bodner-Partom  constitutive  equations  in  STEALTH. 
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FREE  SURFRCE  VELOCITY 


Figure  6.  Comparison  of  Free  Surface  Velocity  Histories 
for  the  Indirect  Approach.  (The  plots  from 
STEALTH  with  and  without  using  B-P  Model  are 
almost  identical.) 


21 


REFERENCES 


1.  Zukas,  J.  "Numerical  Simulation  of  Impact  Phenomena," 

Impact  Dynamics,  Chapter  10,  John  Wiley,  New  York, 
pp. 367-417,  (1981). 

2.  Bodner,S.  R.,  and  Partom,  Y.,  "Constitutive  Equations  for 
Elastic- Viscoplastic  Strain  Hardening  Materials,"  J . 

Applied  Mechanics,  Vol.  *12,  pp.  395-389,  (1  975). 

3.  Walker,  K.  P.,  "Representation  of  Hastelloy-X  Behavior  at 
Elevated  Temperature  with  a  Functional  Theory  of 
Viscoplasticity,"  ASME,  J.  of  Eng.  Mat.  &  Technology,  1981. 

4.  Chaboche,  J.  L.,  "Thermodynamic  and  Phenomenological 
Description  of  Cyclic  Viscoplasticity  with  Damage,"  Trans, 
of  Publicatiuon  No.  1978-3  of  the  Office  National  d’Etudes 
et  de  Researches  Aerospatiales ,  France,  European  Space 
Agency  Technical  Translation  Service,  Publication  No.  ESA- 
548,  May  1979. 

5.  Miller,  A.  K.,  "An  Inelastic  Constitutive  Model  for 
Monotonic,  Cyclic  and  C-Deformat ion :  Part  I  -  Equations 
Development  and  Analytical  Procedures  and  Part  II  - 
Application  to  Type  304  Stainless  Steel,"  ASME,  J .  Eng . 

Mat .  and  Tech . ,  Vol.  98,  pp .  97-113,  (1976). 

6.  Krieg,  R.  D.,  Swearengen,  J.  C. ,  and  Rohde,  R.  W. ,  "A 
Physically  Based  Internal  Variable  Model  for  Rate-Dependent 
Plasticity,"  Inelastic  Behavior  of  Pressure  Vessel  and 
Piping  Components,  ASME,  PVP-PB-028,  pp .  15-28,  (1978). 

7.  "'STEALTH',  Lagrange  Explicit  Finite  Difference  Code  for 
Solids,  Structural  and  Thermohydraulic  Analysis,"  EPRI  NP- 
2080,  November  1981. 

8.  Rajendran,  A.  M.,  and  Grove,  D.  J.,  "Simulations  of  One- 
Dimensional  Plate  Using  STEALTH  -  Finite  Difference  Code," 
UDR-TR-84-1 09 ,  September  1984. 

9.  Bodner,  S.  R.,  and  Merzer,  A.  M.,  "Viscoplastic 
Constitutive  Equations  for  Copper  with  Strain  Rate  History 
and  Temperature  Effects,"  J.  of  Eng.  Mat.  Tech.,  Vol.  100, 
pp.  388-394,  (1978). 

10.  Rajendran,  A.  M. ,  and  Bless,  S.  J.,  "High  Strain  Rate 
Material  Behavior,"  AFWAL-TR-85-4009 ,  Materials  Lab., 
Wright-Patterson  Air  Force  Base,  Ohio  45433-6533,  December 
1985. 


22 


APPENDIX  A 
STEALTH  ADDITIONS 


The  Bodner-Partom  plastic  flow  model  has  been  successfully 
incorporated  into  the  STEALTH  ID  computer  code.  This  task 
required  new  common  block  variables,  additional  input  parameters, 
additions  to  twelve  STEALTH  routines,  and  six  new  subroutines. 

A  total  of  62  new  variables  were  added  to  four  STEALTH 
common  blocks  (/MATARY/,  /MATVAR/,  /ZONARY/,  and  /ZONVAR/) .  The 
new  variables  in  /MATARY/  and  /MATVAR/  are  used  to  store  values 
for  the  input  material-dependent  Bodner-Partom  parameters,  while 
those  in  /ZONARY/  and  /ZONVAR/  are  used  to  store  the  additional 
zone  variables  used  in  the  Bodner-Partom  solution.  Tables  A-1  and 
A-2  describe  these  new  common  variables  in  more  detail.  Five  of 
the  new  zone  variables  (PWK,  EP1 ,  PS1 ,  TS1 ,  and  EFS)  have  been 
given  identifier  numbers  so  that  they  may  be  printed  and/or 
plotted.  However,  the  three  new  zone  variables  relating  to 
plastic  flow  (PWK,  EP1 ,  and  PS1 )  are  only  computed  by  the  Bodner- 
Partom  model,  so  printing  or  plotting  these  variables  is 
meaningless  during  a  "standard”  STEALTH  run  without  using  the  B-P 
model.  Appendix  B  provides  the  identifier  numbers  for  these  new 
zone  variables. 

Five  new  input  record  types  have  been  added  to  the  material 
specification  phase.  Record  type  116  is  used  to  indicate  the 
plastic  flow  model  type  (either  standard  STEALTH  or  STEALTH  with 
the  Bodner-Partom  model),  while  record  types  1 8 1 ,  182,  183,  and 
184  specify  constants,  constraints,  and  convergence  parameters 
required  by  the  Bodner-Partom  plastic  flow  model.  Omission  of 
these  new  input  records  will  result  in  a  "standard"  STEALTH  run. 
However,  if  record  type  116  is  used  to  select  the  Bodner-Partom 
plastic  flow  model  for  a  particular  material,  then  record  types 
1 8 1 ,  182,  and  183  must  also  be  included  for  that  same  material. 

The  convergence  parameters  on  record  type  184  have  default  values, 
so  this  record  is  always  optional.  Appendix  C  describes  the  new 
input  record  types. 
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Two  initialization  routines  (CMNMAT  and  CMNZON)  were  changed 
in  order  to  initialize  all  the  new  common  block  variables. 

INIPTR,  another  initialization  routine,  was  modified  so  that  five 
of  the  new  zone  variables  (PWK,  EP1 ,  PS1 ,  TS1 ,  and  EFS)  could  be 
printed  and  plotted.  Code  was  added  to  three  input  routines 
(MATINP,  MATCHK,  and  MATPRT)  so  that  the  Bodner-Partom  parameters 
could  be  read  in,  checked,  and  printed  out  during  the  material 
specification  phase.  In  addition,  MATPRM  was  changed  so  that  the 
values  of  the  new  variables  in  /MATVAR/  would  always  correspond  to 
the  current  material.  Code  was  added  to  subroutines  ZONOLD  and 
ZONNEW  in  order  to  retrieve  old  values  and  store  new  values  for 
the  new  zone  variables  that  were  added  to  /ZONARY/.  In  subroutine 
ZONSDV,  code  was  added  to  determine  when  the  Bodner-Partom  model 
should  be  used  for  plastic  flow  calculations.  For  plotting 
purposes  only,  subroutines  ZONSTN  and  ZONSTR  were  modified  to 
compute  total  xx  strain  and  effective  stress,  respectively. 
Finally,  six  new  subroutines  were  added  to  the  STEALTH  ID  code  in 
order  to  compute  the  Bodner-Partom  plastic  flow  solution:  BODPAR, 
EPIJ,  SDOTIJ,  SIJDIF,  SI JEST ,  and  SIJNEW.  Table  A-3  summarizes 
the  changes  in  each  STEALTH  routine,  and  Table  A-4  briefly 
describes  each  new  subroutine.  Appendix  D  provides  flowcharts  for 
the  new  routines  that  were  added  to  STEALTH. 
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TABLE  A- 1 


Common  Block 


/MATARY/ 


/MATVAR/ 


/ZONARY/ 


/ZONVAR/ 


MODIFIED  STEALTH  COMMON  BLOCKS 
FOR  BODNER-PARTOM  PLASTIC  FLOW  MODEL 


Additions 


Added  16  array  variables:  MAPLS(IO),  AD0(10), 
ASRN(IO),  AHRMO (10),  AHRMI(IO),  AHRA(IO), 

AZO (10),  AZI(IO),  ADEMAX(IO),  AEFSMN(IO), 
AD2PMN (10),  ATOLSD(IO),  ARLXSD(IO), 
MAITSD(IO),  MACTSD(IO),  ATCFSD(IO).  See 
Table  A-2  for  variable  definitions. 

Added  16  variables:  MPLS,  DO,  SRN,  HRMO, 

HRM1 ,  HRA,  ZO,  Z1 ,  DEMAX,  EFSMIN ,  D2PMIN, 
TOLRSD,  RLAXSD ,  MITRSD,  MCUTSD,  TCFSD.  See 
Table  A-2  for  variable  definitions. 

Added  8  array  variables:  PWK(202),  EP1(202), 
PS1 (202) ,  TS1(202),  EFS(202),  SXXDT(202), 
SYYDT ( 202 ) ,  SZZDT(202).  See  Table  A-2  for 
variable  definitions. 

Added  22  variables:  PWKO,  PWKN,  EP10,  EP1N, 
PS10,  PS1N,  TS1 0 ,  TS1N,  EFSO,  EFSN ,  SXXDTO, 
SXXDTN ,  SYYDTO,  SYYDTN ,  SZZDTO,  SZZDTN, 

SXYDTO ,  SXYDTN ,  SYZDTO,  SYZDTN ,  SXZDTO, 
SXZDTN.  See  Table  A-2  for  variable 
definitions. 
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TABLE  A- 2 


NEW  COMMON  VARIABLES  IN  STEALTH  ID 


VARIABLE 

NAME 

COMMON  BLOCK 
ADDED  TO 

DESCRIPTION 

MAPLS(L) 

/MATARY/ 

Type  of  plastic  flow  model 
for  material  number  L. 

ADO ( L) 

it 

Limiting  value  of  strain  rate 
for  Bodner-Partom  constants 
associated  with  material 
number  L. 

ASRN(L) 

it 

Strain  rate  exponent  for 
material  number  L. 

AHRMO(L) 

IT 

Hardening  rate  parameter  for 
material  number  L. 

AHRM1 (L) 

If 

Hardening  rate  parameter  for 
material  number  L. 

AHRA(L) 

II 

Hardening  rate  parameter  for 
material  number  L. 

AZO(L) 

II 

Initial  value  of  hardness  for 
material  number  L. 

AZ1 (L) 

II 

Maximum  value  of  hardness  for 
material  number  L. 

ADEMAX(L) 

II 

Maximum  allowable  change  in 
the  total  strain  during  one 
time  step  for  material  number 
L. 

AEFSMN(L) 

II 

Minimum  value  of  the  effective 
stress  for  plastic  flow  in 
material  number  L. 

AD2PMN ( L) 

II 

Minimum  value  of  the  second 
invariant  of  the  deviatoric 
plastic  strain  rate  for 
plastic  flow  in  material 
number  L. 

ATOLSD(L) 

It 

Relative  convergence  tol- 

erance  for  the  deviatoric 
stress  computations  for 
material  number  L. 
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ARLXSD(L) 


/MATARY/ 


MAITSD(L)  " 

MACTSD(L)  " 

ATCFSD(L)  " 

MPLS  /MATVAR/ 

DO  " 

SRN  " 

HRMO  " 

HRM1  " 

HR  A  " 

ZO  " 

Z1  " 


Relaxation  factor  used  to 
accelerate  the  convergence  of 
the  deviatoric  stresses  for 
material  number  L. 

Maximum  number  of  iterations 
allowed  during  computation  of 
the  deviatoric  stresses  for 
material  number  L. 

Maximum  number  of  time  step 
cuts  allowed  during  compu¬ 
tation  of  the  deviatoric 
stresses  for  material  number 
L. 

Time  step  cut  factor  used  to 
decrease  the  time  step  length 
when  the  deviatoric  stresses 
fail  to  converge  for  material 
number  L. 

Plastic  flow  model  type. 

Limiting  value  of  strain  rate 
for  Bodner-Partom  constants. 

Strain  rate  exponent. 

Hardening  rate  parameter. 

Hardening  rate  parameter. 

Hardening  rate  parameter. 

Initial  value  of  hardness. 

Maximum  value  of  hardness. 


DEMAX  "  Maximum  allowable  change  in 

the  total  strain  during  one 
time  step. 

EFSMIN  "  Minimum  value  of  the  effective 

stress  for  plastic  flow. 

D2PMIN  "  Minimum  value  of  the  second 

invariant  of  the  deviatoric 
plastic  strain  rate  for 

plastic  flow. 

TOLRSD  "  Relative  convergence  tolerance 

for  the  deviatoric  stress 
computations. 
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RLAXSD 


/MATVAR/ 


MITRSD 

MCUTSD 

TCFSD 

PWK( i) 

EP1 ( i) 

PS1 (  i) 

TS1 ( i) 

EFS( i) 

SXXDT ( i ) 

SYYDT ( i ) 

SZ  ZDT ( i) 

PWKO 
PWKN 
EP 1  0 

EP1N 

PS10 
PS1  N 


it 

ti 

it 

/ZONARY/ 

it 

it 

it 

ti 

ti 

it 

ti 

/ZONVAR/ 

it 

it 

ti 

ti 

ti 


Relaxation  factor  used  to 
accelerate  the  convergence  of 
the  deviatoric  stresses. 

Maximum  number  of  iterations 
allowed  during  computation  of 
the  deviatoric  stresses. 

Maximum  number  of  time  step 
cuts  allowed  during 
computation  of  the  deviatoric 
stresses . 

Time  step  cut  factor  used  to 
decrease  the  time  step  length 
when  the  deviatoric  stresses 
fail  to  converge. 

Total  plastic  work  done  in 
zone  number  i. 

Plastic  xx  strain  rate  in  zone 
number  i . 

Plastic  xx  strain  in  zone 
number  i. 

Total  xx  strain  in  zone 
number  i. 

Effective  stress  in  zone 
nu  mb  e  r  i  . 

Deviatoric  xx  normal  stress 
gradient  in  zone  number  i. 

Deviatoric  yy  normal  stress 
gradient  in  zone  number  i. 

Deviatoric  zz  normal  stress 
gradient  in  zone  number  i. 

Plastic  work  at  old  time. 


Plastic  work  at  new  time. 

Plastic  xx  strain  rate  at 
old  time. 

Plastic  xx  strain  rate  at 
new  time. 

Plastic  xx  strain  at  old  time. 
Plastic  xx  strain  at  new  time. 
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TS1  0 

/ZONVAR/ 

Total  xx  strain  at  old  time. 

TS1N 

tt 

Total  xx  strain  at  new  time. 

EFSO 

ti 

Effective  stress  at  old  time 

EFSN 

ti 

Effective  stress  at  new  time 

SXXDTO 

ti 

Deviatoric  xx  normal  stress 
gradient  at  old  time. 

SXXDTN 

n 

Deviatoric  xx  normal  stress 
gradient  at  new  time. 

SYYDTO 

n 

Deviatoric  yy  normal  stress 
gradient  at  old  time. 

SYYDTN 

n 

Deviatoric  yy  normal  stress 
gradient  at  new  time. 

SZZDTO 

n 

Deviatoric  zz  normal  stress 
gradient  at  old  time. 

SZZDTN 

ii 

Deviatoric  zz  normal  stress 
gradient  at  new  time. 

SXYDTO 

ii 

Deviatoric  xy  shear  stress 
gradient  at  old  time. 

SXYDTN 

ii 

Deviatoric  xy  shear  stress 
gradient  at  new  time. 

SYZDTO 

ii 

Deviatoric  yz  shear  stress 
gradient  at  old  time. 

SYZDTN 

ii 

Deviatoric  yz  shear  stress 
gradient  at  new  time. 

SXZDTO 

ii 

Deviatoric  xz  shear  stress 
gradient  at  old  time. 

SXZDTN 

Deviatoric  xz  shear  stress 
gradient  at  new  time. 
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TABLE  A-3 


ADDITIONS  TO  EXISTING  STEALTH  ROUTINES 


Routine 


Additions 


CMNMAT 

CMNZON 

INIPTR 


MATCHK 

MATINP 

MATPRM 

MATPRT 

ZONNEW 

ZONOLD 

ZONSDV 

ZONSTN 

ZONSTR 


Added  code  to  initialize  the  new  variables  that 
were  added  to  common  blocks  /MATARY/  and  /MATVAR/. 

Added  code  to  initialize  the  new  variables  that 
were  added  to  common  blocks  /ZONARY/  and  /ZONVAR/. 

Added  5  new  variable  names  (PWK,  EP1  ,  PS1  ,  TS1  ,  and 
EFS)  to  the  NAMVAR  array  data  statement  (at 
locations  77,  78,  79,  80,  and  90  in  the  array). 

Also  added  code  to  store  addresses  and  data  types 
of  these  variables  (in  IDPTR  and  IDATT  arrays)  so 
that  they  can  be  printed  and/or  plotted. 

Added  code  to  check  Bodner-Partom  plastic  flow 
model  input  parameters  (from  record  types  181,  182, 
183,  and  184)  for  consistency  and  validity. 

Added  code  to  read  5  new  types  of  MAT  input  records 
(116,  181,  182,  183,  and  184).  The  116  input 
record  specifies  the  plastic  flow  model  type,  while 
the  4  other  input  records  define  the  parameters 
required  by  the  Bodner-Partom  plastic  flow  model. 

Added  code  to  define  the  new  variables  that  were 
added  to  common  /MATVAR/.  These  variables  contain 
the  Bodner-Partom  plastic  flow  model  parameters  for 
the  current  material. 

Added  code  to  echo  the  Bodner-Partom  plastic  flow 
model  input  parameters  (when  appropriate)  for 
each  material. 

Added  code  to  store  new  values  of  the  zone 
variables  related  to  the  Bodner-Partom  plastic  flow 
model.  These  values  are  stored  in  the  new  array 
variables  that  were  added  to  common  /ZONARY/. 

Added  code  to  define  old  values  of  the  zone 
variables  related  to  the  Bodner-Partom  plastic 
flow  model.  These  values  are  retrieved  from  the 
new  array  variables  that  were  added  to  common 
/ZONARY/. 

Added  code  to  compute  deviatoric  stress  gradients 
for  elastic  flow  and  to  call  Bodner-Partom  routines 
for  plastic  flow  (when  appropriate). 

Added  code  to  compute  new  total  xx  strain. 

Added  code  to  compute  new  effective  stress. 
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TABLE  k-H 


Routine 

BODPAR 

EPIJ 

SDOTIJ 

SIJDIF 

SI JEST 

SIJNEW 


NEW  STEALTH  ROUTINES 
FOR  BODNER-PARTOM  PLASTIC  FLOW  MODEL 

_ Description _ 


Controls  the  Bodner-Partom  plastic  flow  model 
computations.  This  routine  divides  the  STEALTH 
time  step  into  smaller  steps  to  facilitate  con¬ 
vergence  of  the  Bodner-Partom  calculations.  Output 
from  this  routine  consists  of  the  deviatoric 
stresses  and  stress  gradients,  the  plastic  work 
done  (so  far),  the  plastic  xx  strain  rate,  and  the 
plastic  xx  strain. 

Uses  the  Bodner-Partom  viscoplastic  state  variable 
theory  to  compute  the  plastic  strain  rates  and  the 
second  invariant  of  the  deviatoric  plastic  strain 
rate . 

Computes  the  deviatoric  stress  gradients. 

Computes  the  maximum  relative  difference  between 
estimated  and  computed  values  of  the  deviatoric 
stresses. 

Computes  initial  estimates  for  the  new  deviatoric 
stresses. 

Computes  the  new  deviatoric  stresses. 
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APPENDIX  B 


STEALTH  VARIABLE  IDENTIFIERS  (IDS)* 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

TIM 

TIE 

IDE 

THE 

TEE 

TEG 

TXM 

TYM 

TXK 

TYR 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

XPN 

XVL 

XAC 

YPN 

YVL 

YAC 

ZPN 

ZVL 

ZAC 

DST 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

GXK 

GYK 

GZK 

XHF 

YHP 

ZHF 

TZM 

TZK 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

NOR 

NBM 

NBT 

NBV 

ACT 

MPN 

IND 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

DLL 

RLV 

COM 

QDA 

VSR 

DNS 

ZMS 

TRV 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

ZIZ 

ZDE 

ZHE 

ZSE 

ZKE 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

TMP 

CON 

SHC 

PRH 

AVS 

SSP 

BFS 

IGN 

71 

72 

73 

74 

75 

i 

76  ; 

77 

78 

79 

80 

TXX 

TYY 

TZZ 

TXY 

TYZ 

TXZ  ! 

PWK 

EP1 

PS1 

TS1 

81 

82 

83 

84 

85 

86 

87 

88 

89 

!  90 ! 

SXX 

SYY 

SZZ 

SXY 

SYZ 

SZX 

YLD 

SHR 

jEFSj 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

CXI 

EX2 

EX3 

EX4 

EX5 

EX6 

The  Identifiers  ere  index  (address)  values  for  the  arrays  IDPTR  (variable 
pointer)  and  IDATT  (variable  attribute)  found  in  subroutine  INTPTR. 
Variables  1-10  are  global  variables,  while  the  rest  are  either  grid-point 
or  zone-related.  Variable  definitions  may  be  found  in  Appendix  "Conven¬ 
tions"  of  Volume  3,  "Programmer's  Manual." 
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APPENDIX  C 


ADDITIONAL  STEALTH  INPUT  RECORDS 


PLASTIC  FLOW  MODEL  TYPE 


The  116  input  record  specifies  the  type  of 
plastic  flow  model  for  material  MPN.  If  the 
116  input  record  is  omitted,  the  default 
value  is  used.  The  116  input  record  may  be 
used  at  restart  time  to  change  the  plastic 
flow  model  type. 


Columns  11  21  31  Hi  51  61  71  80 

JmPN  I  MPLS  J  —  |  —  J  —  J  —  J - J 

Fields  1  2  3  H  5  6  7 

Input 

Field  Variable 

Number  Name  _ Description  of  Input  Variable _ _ 

1  MPN  Material  property  number.  See  Field 

Number  1  on  TYPE  111  record. 

2  MPLS  Plastic  flow  model  type.  The  value  for 

MPLS  identifies  the  type  of  plastic  flow 
model  for  material  MPN.  The  value  of  MPLS 
must  be  integral  and  must  be  greater  than 
or  equal  to  1  and  less  than  or  equal  to  2. 
If  MPLS  is  out  of  this  range,  a  RGEERR 
message  is  written  on  the  print  and  the 
message  files  and  the  program  terminates 
immediately  in  subroutine  MATINP.  (See 
Appendix  "Messages.")  The  following 
definitions  apply: 

=  1.0  standard  STEALTH  (default) 

=2.0  Bodner-Partom 

Parameters  required  by  the  Bodner-Partom 
model  are  entered  on  input  records  1 8 1 , 
182,  183,  and  1 84 .  There  are  no 
additional  parameters  associated  with  the 
standard  STEALTH  model. 

3~7  -  Fields  3~7  must  be  left  blank  or  the 

FLDCHK  message  is  written  on  the  print  and 
the  message  files  and  the  program 
terminates  at  the  end  of  the  GENERATOR 
phase  group.  (See  Appendix  "Messages.") 
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BODNER-PARTOM  STRAIN  RATE  PARAMETERS 


page  1  of  1 


TYPE 

181 


The  1 81  input  record  specifies  the  strain 
rate  parameters  (for  material  MPN)  required 
by  the  Bodner-Partom  plastic  flow  model.  The 
1 81  input  record  is  absolutely  required  when 
the  Bodner-Partom  model  is  specified  (MPLS=2 
on  the  1 1 6  input  record);  it  is  ignored  when 
the  standard  STEALTH  plastic  flow  model  is 
specified  (MPLS=1  on  the  116  input  record). 

The  1 8 1  record  may  be  used  at  restart  time  to 
change  the  Bodner-Partom  strain  rate  parameters. 


START 

RESTART 


0 

MAXMAT 

0 

MAXMAT 

Columns 


11  21  31  Ml  51  61  71  80 


MPN 

DO 

SRN 

— 

— 

— 

— 

Fields 

1 

Field 

Number 

Input 

Variable 

Name 

1 

MPN 

2 

DO 

3 

SRN 

4-7 

3  4  5  6  7 


Description  of  Input  Variable 


Material  property  number.  See  Field 
Number  1  on  TYPE  111  record. 

Limiting  value  of  strain  rate.  The 
value  of  DO  must  be  greater  than  0.0. 

Strain  rate  exponent.  The  value  of 
SRN  must  be  greater  than  0.0. 

Fields  4-7  must  be  left  blank  or  the 
FLDCHK  message  is  written  on  the  print 
and  the  message  files  and  the  program 
terminates  at  the  end  of  the  GENERATOR 
phase  group.  (See  Appendix  "Messages.") 
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BODNER-PARTOM  HARDENING  PARAMETERS 


page  1  of  1 


TYPE 

182 


The  182  input  record  specifies  the  material 
hardening  parameters  (for  material  MPN)  required 
by  the  Bodner-Partom  plastic  flow  model.  The 
182  input  record  is  absolutely  required  when 
the  Bodner-Partom  model  is  specified  (MPLS=2 
on  the  116  input  record);  it  is  ignored  when 
the  standard  STEALTH  plastic  flow  model  is 
specified  (MPLS=1  on  the  116  input  record). 

The  182  record  may  be  used  at  restart  time 
to  change  the  Bodner-Partom  material 
hardening  parameters. 


START 

RESTART 


0 

MAXMAT 

0 

MAXMAT 

Columns 


Fields 


11  21  31  41  51  61  71  80 


|"mpn 

HRMO 

HRM1  | 

HRA 

zo 

Z1 

1 

2 

3 

H 

5 

6 

7 

Input 

Field  Variable 

Number  Name  _ Description  of  Input  Variable 


1  MPN  Material  property  number.  See  Field 

Number  1  on  TYPE  111  record. 

2  HRMO  Hardening  rate  parameter,  m  .  The  value 

of  HRMO  must  be  greater  than  0.0. 


3  HRM1  Hardening  rate  parameter,  m  .  The  value 

of  HRM1  must  be  greater  than  or  equal  to 

0.0. 

4  HRA  Hardening  rate  parameter,  a.  The  value  of 

HRA  must  be  greater  than  or  equal  to  0.0. 


5  Z0 


Initial  value  of  hardness.  The  value  of 
Z0  must  be  greater  than  0.0. 


6  Z1  Maximum  value  of  hardness.  The  value  of 

Z1  must  be  greater  than  0.0. 

7  —  Field  7  must  be  left  blank  or  the  FLDCHK 

message  is  written  on  the  print  and  the 
message  files  and  the  program  terminates 
at  the  end  of  the  GENERATOR  phase  group. 
(See  Appendix  "Messages.") 
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BODNER-PARTOM  CONSTRAINTS 


page  1  of  1 


TYPE 

183 


The  183  input  record  specifies  constraints 
(for  material  MPN)  required  by  the  Bodner-Partom 
plastic  flow  model.  The  183  input  record  is 
absolutely  required  when  the  Bodner-Partom 
model  is  specified  (MPLS=2  on  the  116  input 
record);  it  is  ignored  when  the  standard  STEALTH 
plastic  flow  model  is  specified  (MPLS=1  on  the 
116  input  record).  The  183  record  may  be  used 
at  restart  time  to  change  the  Bodner-Partom 
constraints. 


START 

RESTART 


0 

MAXMAT 

0 

MAXMAT 

Columns 

11 

21 

31  41  51 

61  71 

|  MPN 

|  DEMAX 

|  EFSMIN  |  D2PMIN  |  — 

- 1  -  1  - 

Fields 

1 

2 

3  4  5 

6  7 

Field 

Number 

Input 

Variable 

Name 

Description  of 

Input  Variable 

1  MPN  Material  property  number.  See  Field 

Number  1  on  TYPE  111  record. 


2  DEMAX 


3  EFSMIN 


4  D2PMIN 


5-7 


Maximum  allowable  change  in  the  total 
strain  during  one  time  step.  The  value 
of  DEMAX  must  be  greater  than  0.0. 

Minimum  value  of  the  effective  stress  for 
plastic  flow.  The  value  of  EFSMIN  must  be 
greater  than  0.0. 

Minimum  value  of  the  second  invariant 
of  the  deviatoric  plastic  strain  rate  for 
plastic  flow.  The  value  of  D2PMIN  must  be 
greater  than  0.0. 

Fields  5-7  must  be  left  blank  or  the 
FLDCHK  message  is  written  on  the  print  and 
the  message  files  and  the  program 
terminates  at  the  end  of  the  GENERATOR 
phase  group.  (See  Appendix  "Messages.") 
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BODNER-PARTOM  CONVERGENCE  PARAMETERS 


page  1  of  2 


TYPE 

184 


The  184  input  record  specifies  convergence 
parameters  (for  material  MPN)  required  by 
the  Bodner-Partom  plastic  flow  model. 

Default  values  are  used  if  the  184  input 
record  is  omitted  when  the  Bodner-Partom 
model  is  specified  (MPLS=2  on  the  116  input 
record) .  This  record  is  ignored  when  the 
standard  STEALTH  plastic  flow  model  is  specified 
(MPLS=1  on  the  116  input  record).  The  184 
input  record  may  be  used  at  restart  time  to 
change  the  Bodner-Partom  convergence  parameters. 


START 

RESTART 


0 

MAXMAT 

0 

MAXMAT 

Columns 

11  21 

31  41 

51  61  71 

|  MPN  | TOLRSD 

|  RLAXSD  J  MITRSD 

MCUTSD  |  TCFSD  |  — 

Fields 

1  2 

3  4 

5  6  7 

Field 

Number 

Input 

Variable 

Name 

Description 

of  Input  Variable 

1 

MPN 

Material  property  number.  See  Field 

Number  1  on  TYPE  111  record. 


2  TOLRSD 


3  RLAXSD 


4  MITRSD 


Relative  convergence  tolerance  for  the 
deviatoric  stresses.  The  value  of  TOLRSD 
must  be  greater  than  0.0.  If  this  field 
is_^eft  blank,  the  default  value  of  1 .0  x 
10  is  assumed. 

Relaxation  factor  used  to  accelerate  the 
convergence  of  the  deviatoric  stresses. 
The  value  of  RLAXSD  must  be  greater  than 
0.0  and  less  than  or  equal  to  1.0.  If 
this  field  is  left  blank,  the  default 
value  of  1.0  is  assumed. 

Maximum  number  of  iterations  allowed 
during  computation  of  the  deviatoric 
stresses.  The  value  of  MITRSD  must  be 
integral  and  greater  than  or  equal  to  1 . 
If  this  field  is  left  blank,  the  default 
value  of  5  is  assumed. 
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BODNER-PARTOM  CONVERGENCE  PARAMETERS 


page  2  of  2 


Columns 

Fields 

Field 

Number 

5 

6 


TYPE 

184 


11  21  31  41  51  61  71  80 


|  MPN 

|tolrsd 

|  RLAXSD 

MITRSD 

MCUTSD 

TCFSD  j 

EJ 

1 

2 

3 

4 

5 

6 

7 

Input 

Variable 

Name  _ Description  of  Input  Variable 


MCUTSD  Maximum  number  of  time  step  cuts  allowed 

during  computation  of  the  deviatoric 
stresses.  The  value  of  MCUTSD  must  be 
integral  and  greater  than  or  equal  to  1 . 
If  this  field  is  left  blank,  the  default 
value  of  10  is  assumed. 


TCFSD  Time  step  cut  factor  used  when  the 

deviatoric  stresses  fail  to  converge.  The 
value  of  TCFSD  must  be  greater  than  1.0. 
When  a  time  step  cut  is  made,  the  length 
of  the  new  time  step  is  equal  to  the  old 
time  step  divided  by  TCFSD;  i.e.,  TCFSD  = 
2.0  causes  the  time  step  length  to  be  cut 
in  half  at  each  time  step  cut.  If  this 
field  is  left  blank,  the  default  value  of 
2.0  is  assumed. 


Field  7  must  be  left  blank  or  the  FLDCHK 
message  is  written  on  the  print  and  the 
message  files  and  the  program  terminates 
at  the  end  of  the  GENERATOR  phase  group. 
(See  Appendix  "Messages.") 
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ZONSDV 


COMPUTES 

DEVIATORIC 

STRESSES 


APPENDIX  D 

FLOWCHARTS  ON  NEW  SUBROUTINES 


EP1J 


USES  lODrei-PARTOM 
theory  to  compute 

PLASTIC  STRAIN 
RATES 


B0DPA1 


CONTROLS 
CALCULATION  OR 
OEVUTORIC 
STRESSES  DURING 
PLASTIC  PLOW 


SI JEST 


COMPUTES  INITIAL 
ESTIMATES  TOR  NEW 
DEVIATOR I C  STRESSES 


SCOT I J 


COMPUTES 

DEVIATORIC 

STRESS 

gradients 


COMPUTES  NEW 

DEVIATORIC 

STRESSES 


EPIJ 

SDOTIJ 

COMPUTES 

DEVIATORIC 

STRESS 

GRADIENTS 


USES  BODNER-PARTOM 
THEORY  TO  COMPUTE 
PLASTIC  STRAIN  RATES 


L 

cr  mrv 

a  UB»» 

COMPUTES  MAXIMUM 
RELATIVE  DIFFERENCE 
BETWEEN  ESTIMATED 
AND  COMPUTED  VALUES 
OF  DEVIATORIC  STRESSES 


EPCJ 


USeS  BODNER-PARTOM 
THEORY  to  COMPUTE 
PLASTIC  STRAIN  RATES 


EPU 


USES  BODNER-PARTOM 
THEORY  TO  COMPUTE 
PLASTIC  STRAW  RATES 


FLOW  DIAGRAM  FOR  BODNER-PARTOM  PLASTIC 
FLOW  MODEL  IN  STEALTH  ID 
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Vi*  cowBon  block/ZONVAR/ 


BODPAR  Flowchart 


40 


INDICATE  SUCCESSFUL 
CONVERGENCE  BY 

SETTING  CONVERGENCE 

FLAG  EQUAL  TO  1 

f  START  V- - - 

ZET 


(SIJNEW  requires  Cj,  C2.  C3,  P* 
At,,  UP.  S°]d,  s2jd .  and  the 
I  input  convergence  parameters 


CALL  SIJEST 


I  SI  JEST  computes  the  lntial 
Idevlatoric  stress  estimates. 


r 


LOOP  OVER 
MAXIMUM  NUMBER 
OF  ITERATIONS 
ALLOWED 


CALL  SDOTIJ 


1 


b— £ 


SDOTIJ  computes  the  new 
deviator ic  stress  gradients.  S 


ij 


COMPUTE  NEW 
DEVIATORIC 

_ j 

STRESSES 

1 

.new  eold  ‘new  A 
Sij  'SiJ  +Sij  4t. 


1 - £ 


SIJDIF  computes  the  relative 
Indifferences  between  and  sj*w 


I  Is  the  maximum  relative 
— m  J  difference  (from  SIJDIF) 
\  within  the  convergence 
I  tolerance? 


USE  RELAXATION 
FACTOR  TO  DEFINE 
DEVIATORIC  STRESS 
ESTIMATES  FOR  NEXT 
ITERATION 


-{■ 


s"'  •  -c  * 


1 - 


END  OF 

ITERATION  LOOP 


INDICATE  CONVERGENCE 
FAILURE  BY  SETTING 
CONVERGENCE  FLAG 
EQUAL  TO  0 


i 


^  RETURN  -  -  - 


SIJNEW  returns  S. . 

-new  J 

U  ’ 

^lag  via  the  argument  list 


SIJNEW  Flowchart 


SDOTIJ  use*  ch* 

Co  compute  new  devletorlc 
seres*  gradient  estimate*,  S 


eat 

U 


SIJEST  Flowchart 
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SDOTIJ  Flowchart 
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r 

H 

— ,  Dq.  «.  n.  Zo,  Zlf  (a#ff)«lB, 
and  (D^)bIo 


EPIJ  raqulraa  S^,  Wp, 


INITIALIZE 

PLASTIC 

FLOW  PARAMETERS 


r 

_  _____  -  o.o 

"  ~ tv 00 


- ^  °aff  '  2  E  SJj 


YES  <a 

^FlOW  OEflEMTEL^— 
II  ELASTtl 


- T1*  °.f(  « 


COMPUTE  THE  INTERNAL 
STATE  VARIABLE  THAT 
GOVERNS  RESISTANCE 
TO  PLASTIC  FLOW 


I 


- y-h-  <w'*P 


COMPUTE  SECOND 
INVARIANT  OF  THE 
OERIVATORIC  PLASTIC 
STRAIN  RATE 


- {' 


o?  -  o2  .-[(z2/^ff,n 

2  o 


COMPUTE 

PLASTIC 

RATES 

THE 

STRAIN 

b— ^ 


la  0\  <  (Dinin’ 


'<}  *  xsu 


/  *  \  r 

f  arrow  J—  —  —  — 


EPIJ  raturna  dJ  and  e ^ 
via  tha  arguaant  Uat 


EPIJ  Flowchart 
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APPENDIX  E 


INPUT  DATA  FOR  STEALTH  TO 
ANALYZE  A  TENSILE  BAR 


TTL  ID 

THIN  ROD. 

V=10  M/S, 

B-P,  RATE 

INDEPENDENT. 

F'RB 

SYM 

6.0000 

GRD 

1.0000 

0.0000 

1 1 . 0000 

END 

NAT 

111 

1 . 0000 

1 . 0000 

112 

1 . 0000 

2.0000 

2.0000 

116  PL8M0D 

1 . 0000 

2. 0000 

121 

1 . 0000 

7.8560 

122 

1 . 0000 

1 . 0000 

1 . 5400 

123 

1 . 0000 

-10.0000 

10.0000  . 

33663366 

154.00 

132 

1 . 0000 

999 . 0000 

134 

1 . 0000 

0. 7940 

131  BPSTRN 

1 . 0000 

100.0000 

6.5000 

182  BPHARD 

1 . 0000 

200.0000 

0.0000 

0.0000 

0.0150 

133  BPCN3T 

1 . 0000 

1 . OE-4 

1 . 0E-5 

1.0E-10 

134  BPCONV 

1 , 0000 

1.0E-4 

1 . 0000 

10.0000 

10.0000 

END 

OPT 

211 

1 . 0000 

1 . 0000 

212 

1 . 0000 

-5.0000 

0.5000 

2.0000 

212 

2.0000 

50.0000 

0.5000 

3.0000 

221 

1 . 0000 

1 . 0000 

1 1 . 0000 

241 

1 . 0000 

10.0000 

0.0000 

10.0000 

END 

ZON 

311 

1 . 0000 

1 . 0000 

1 1 . 0000 

321 

1 . 0000 

1 . 0000 

322 

1 . 0000 

1 . 0000 

END 

BDY 

421 

1 . 0000 

5.0000 

5.0000 

421 

1 1 . 0000 

4.0000 

5.0000 

1 . 0000 

422 

1 . 0000 

1 . 0000 

441 

1 . 0000 

1 . 0000 

0.0000 

1 . 0E+4 

442 

1 . 0000 

0.0010 

0.0000 

END 

TIN 

511 

0.0100 

512 

0.0050 

10.0000 

521 

1000.00 

1000.00 

END 

EDT 

621  PRTEDT 

1 . 0000 

999.000 

1000.00 

1 . 0000 

622  PRTPTS 

1 . 0000 

1 1 . 0000 

623  PRTIDS 

1 1 . 0000 

12.0000 

45.0000 

71.0000 

72.0000 

671  PLTEDT 

1 . 0000 

0.0000 

1000.000 

1 . 0000 

672  PLOT 

1 . 0000 

1 . 0000 

71.0000 

80.0000 

6.0000 

672 

2.0000 

l . 0000 

71.0000 

1 . 0000 

6.0000 

672 

3.0000 

1 . 0000 

80.0000 

1 . 0000 

6.0000 

END 

END 


1 . 0000 


0.0150 

2.0000 


64.0000  66.0000 
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APPENDIX  F 


INPUT  DATA  FOR  STEALTH  TO 
ANALYZE  A  PLATE  IMPACT  TEST 


TTL 

1 

-D  PLATE 

IMPACT,  RATE 

INDEPENDENT.  IMP 

VEL  *  300 

M/S.  3MM 

PRB 

MIN 

MINUM 

1 . 0E-08 

1.0E-06 

GRD 

FLYER 

1 . 0000 

1 . 0000 

31.0000 

GRD 

TARG 

2.0000 

1 . 0000 

61.0000 

END 

MAT 

111 

EOS 

1 . 0000 

1 . 0000 

112 

STREN 

l . 0000 

2.0000 

2.0000 

2.0000 

116 

PLSMOD 

1 . 0000 

2.0000 

121 

DENS 

1 . 0000 

78.5600 

122 

BULK 

1 . 0000 

1 . 0000 

1540.0000 

t32 

YIELD 

l . 0000 

999.0000 

134 

SHEAR 

1 . 0000 

794.0000 

136 

SPALL 

1 . 0000 

-999 . 0000 

181 

BPSTRN 

1 . 0000 

100.0000 

6.5000 

182 

BPHARD 

1 . 0000 

0.2000 

0. 0000 

0.0000 

15.0000 

15.0000 

183 

BPCNST 

1 . 0000 

1 . OE-4 

1.0E-5 

O 

w 

o 

• 

H 

134 

BPCONV 

1 . 0000 

1 . OE-4 

1 . 0000 

10.0000 

10.0000 

2.0000 

END 

GPT 

211 

I L »  XL 

1 . 0000 

1 . 0000 

-3.0000 

211 

2.0000 

32.0000 

0.0000 

221 

IL,  IR 

1 . 0000 

1 . 0000 

31.0000 

221 

2.0000 

32.0000 

92.0000 

241 

XL,  XR 

1 . 0000 

30.0000 

-3.0000 

0.0000 

241 

2.0000 

60 ■ 0000 

0.0000 

6.0000 

END 

ZON 

311 

IL,  IR 

1 . 0000 

1 . 0000 

31.0000 

311 

2.0000 

32.0000 

92.0000 

321 

MPN 

1 . 0000 

1 . 0000 

321 

2.0000 

1 . 0000 

322 

RLV 

1 . 0000 

1 . 0000 

322 

2.0000 

1 . 0000 

331 

XVL 

1 . 0000 

0.3000 

END 

BOY 

421 

LBCFLR 

1 . 0000 

3.0000 

5.0000 

421 

RBCFLR 

31.0000 

7.0000 

5.0000 

421 

LBCTRG 

32.0000 

7.0000 

5.0000 

421 

RBCTRG 

92.0000 

3 . 0000 

5.0000 

END 

TIM 

511 

TIMSTP 

0.0020 

512 

STPLIM 

0.0019 

10.0000 

513 

STPFCT 

1 . 2000 

0.9000 

521 

TIMLIM 

2.6000 

1000.0000 

END 

EDT 

671 

PLTEDT 

1 . 0000 

0.0000 

2.6000 

0.0100 

672 

PLOT 

1 . 0000 

1 . 0000 

12.0000 

1 . 0000 

92.0000 

672 

2.0000 

1 . 0000 

71.0000 

1 . 0000 

62.0000 

END 

END 


>  6MM 
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